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We show theoretically that elongated vortex-antivortex dipoles can be created controUably in 
trapped Bose-Einstein condensates, using known experimental techniques. Vortex dipoles of suf- 
ficient length are unstable and cascade into slow vortex rings which ultimately decay via sound 
emission. This instability of antiparallel vortex line elements, which self-generates Kelvin waves 
on vortex loops and in trapped atomic gases, may play a role in bridging the Kelvin- wave and 
Kolmogorov-Richardson cascades of quantum turbulence. 
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Contrails in the sky left behind by aircraft can reveal 
the pair of counter-circulating wing-tip vortices gener- 
ated in the wake of the plane [1]. A large proportion 
of the energy required to keep an aircraft airborne is 
consumed in the continuous generation of such wing-tip 
vortices. These powerful eddies create a bottleneck at 
airports due to the hazard they impose on aircrafts fly- 
ing in their vicinity. The Crow instability mechanism 
seeded by atmospheric turbulence is considered to be a 
major agent responsible for breaking up these coherent 
and long-lived wing-tip vortex pairs [21 [3] . In water vor- 
tex dipoles, an additional short-wave instability has been 
observed to grow due to a three-dimensional elliptic in- 
stability mechanism in an antisymmetric mode together 
with the long- wave Crow instability [1] . A quantum ana- 
logue of the Crow instability [51 |B] may occur when bodies 
are dragged through superfluids creating vortices whose 
circulation is quantized [T . 

Bose-Einstein condensates (BECs) are versatile quan- 
tum liquids exhibiting rich superfluid dynamics [8]. 
Quantized vorticity and persistent currents are the hall- 
mark of superfluidity in BECs [rOHTS] . The analog of the 
Kadomtsev-Petviashvili [IB] or 'snake' instability of soli- 
tons has been observed to lead to generation of vortices 
and vortex rings in BECs [T7H22] . In the limit of antipar- 
allel vortices it further transforms into the Crow insta- 
bility [S] . In pancake condensates vortex-antivortex pairs 
or vortex dipoles can be nucleated via active stimulation 
[23j . For suitable parameter regimes such vortex shed- 
ding from a moving obstacle [25] is theoretically predicted 
to exhibit a Benard-von Karman vortex street structure 
[21]. Vortex dipoles can also form spontaneously via 
the Berezinskii-Kosterlitz-Thouless mechanism in finite- 
temperature systems in which entropy may cover the en- 
ergy cost of pair creation ;26ii30i[ . Furthermore, a quench 
through a BEC phase transition facilitates stochastic for- 
mation of vortex dipoles via the Kibble-Zurek scenario as 
observed in recent experiments |31l 132) . 

Solitary vortex dipoles have been created and ob- 
served experimentally in oblate condensates. However, 
the three-dimensional Crow instability is manifestly ab- 
sent in sufficiently fiat condensates where axial degrees 
of freedom play no role in vortex dynamics. Since spon- 
taneous formation of long three-dimensional vortices is 



suppressed, an active pair-creation method is needed to 
isolate and observe three-dimensional instabilities of in- 
teracting vortex lines. Here we propose such a method to 
generate long vortex dipoles in trapped quantum degen- 
erate gases and study their superfluid decay dynamics. 

Bose-Einstein condensates may be represented by a 
macroscopic wavefunction 0(r, t) whose evolution is gov- 
erned by the Gross-Pitaevskii equation 
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where the constant a is the s-wave scattering length of 
particles with mass m, and Vcxt(r) is an external poten- 
tial used to confine and manipulate the atoms. We first 
consider an initial state wave function (p{r, 0) = {x+d/2+ 
e cos(fcz) -f- iy) {x — d/2 — e cos(fcz) — iy) which uses Carte- 
sian coordinates to represent two counter-circulating vor- 
tex lines separated by distance d. The vortices are sinu- 
soidally perturbed with an amplitude e and a wave vector 
fc in a symmetric mode about their equilibrium positions, 
see Fig. 1(a). This vortex dipole generates a mutual in- 
duction field which causes the pair to travel with a speed 
inversely proportional to d in the positive y direction. In 
addition to the translational motion, each perturbed vor- 
tex spins about its own axis. This self-induced spinning 
motion is due to the curvature of the vortex and comes in 
the form of Kelvin waves of growing amplitude 33, 34]. 
In the Crow instability mode the vortices become phase- 
locked at a certain angle when the total induced velocity 
field stops the vortices from spinning. Eventually parts of 
the vortices overlap, forming a chain of vortex reconnec- 
tions, shown in Fig. 1(c), which break the vortex dipole 
into a sequence of vortex loops, see Fig. 1(d). The dy- 
namics illustrated in Fig.l is obtained by considering a 
homogeneous (Vcxt = 0) non-interacting (a = 0) system, 
formally integrating Eq. (fTl) and truncating the power 
series expansion of the resulting matrix exponential to 
second order, yielding an analytically soluble model. Vor- 
tices are visualized in Fig. 1 by plotting isosurfaces of the 
function |0(r,i)p. 

Next we model the full nonlinear dynamics of an elon- 
gated vortex dipole embedded in an inhomogeneous back- 
ground by numerically solving the full Gross-Pitaevskii 
equation. We consider N = J |0(r)pdr = 8 x 10^ Bose- 
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FIG. 1. Evolution of a sinusoidally perturbed vortex-antivortex dipole. (a) The superflow around the vortices is denoted by 
the circulation vectors k,± and the oriented circles on top of the subfigure indicate the rotational motion of vortex line elements 
in the direction opposite to the circulating superflow. The pair translates in the y direction shown by the arrow, due to the 
mutual induction field, (b) The amplitude of the instability grows, (c) Vortex reconnections occur at locations where vortices 
come in contact with each other, (d) The initial vortex pair has broken into a sequence of vortex loops, (e) Side view of frame 
d. The projection of the vortex lines on the z = —2d plane are shown on the bottom of each subfigure a-d which are plotted, 
respectively, for times i/io =0.02, 0.04, 0.05, and 0.06, where the unit of time to = AmcP /h. Parameters described in the text 
are e = d/200 and k — 5/d. 



Einstein condensed ®''Rb atoms in an anisotropic har- 
monic potential 14xt(r) = rnuj^{x'^ 



y 



Xlz^)/2, where 
A2 — cjz/^ = 0.2, and oj = 2tt x lOOHz. The conden- 
sate heaUng length ^0 = 0.15 fj,m characterizes the size of 
the vortex core while for our parameters the harmonic 
oscillator length oq = ^/h/muj = 6.6^0 ■ We prepare the 
initial state by phase imprinting two straight e = vor- 
tex lines of opposite circulation in the condensate wave- 
function. The cores of the vortices are placed at loca- 
tions {xo = ±d/2,j/o = — 2ao) and are oriented along the 
z axis. The ends of the dipole are deliberately short- 
circuited at zo = ±20ao, topologically forming a highly 
anisotropic vortex ring, shown in Fig. 2(a). Note that 
vortices described by a single analytical complex func- 
tion cannot have free ends and will always form closed 
loops, although in practice quantum and/or thermal fluc- 
tuations may provide an effective boundary for the con- 
densate where vortices may terminate. 

The mother vortex ring is propelled forward in the y 
direction by the self-induced superflow. The inhomoge- 
neous condensate density causes the vortex line elements 
in the lower density regions to travel faster than those in 
higher density regions. Due to this differential velocity 
the ends of the dipole bend, self-generating Kelvin waves 
which then propagate along the vortex lines symmetri- 
cally from both ends toward the center of the conden- 
sate, see Fig. 2(a) and Fig. 2(b). These Kelvin waves 
catalyse the growth of instabilities on the vortex dipole. 
A train of daughter vortex loops is created when the 
unstable modes have grown in amplitude to close the 
gap separating the vortices, igniting reconnections at the 
nodes of the most unstable modes, shown in Fig. 2(b) 
and Fig. 2(c). This violent process leaves the daughter 
vortex loops ringing (carrying Kelvin wave excitations) 
transforming them into slow vortex rings |35l 136] . These 
daughter vortex loops subsequently undergo several re- 



TABLE I. Instability parameters for different initial vortex 
dipole separations d. Mean distance dr separating the vortex 
lines at time r is shown in parenthesis. 



rf(€o) 
26(11) 
20(10) 
12(7) 
9(6) 
5(4) 



Ar(go) 
20 
14 
12 
12 
12 



55° 
43° 
40° 
22° 
18° 



r(ms) 
12^ 
9 
6 
4 
2 



combinations and reconnections and ultimately convert 
their energy into sound waves. If the size of the gener- 
ated vortex loops is large enough, they will penetrate the 
condensate surface. This can result in rotation of the axis 
of vorticity by 90 degrees, from initially being oriented 
along the z axis to lying along the x axis, cf. Fig. 2(a) 
and Fig. 2(d). Table I summarizes the dependence of 
the instability on the separation d{^o) between the vor- 
tices. The wavelength of the fastest growing mode At-(^o) 
and the angles /3r and ar (cf. Fig. 1(d) and Fig. 1(e)), 
which in our case are approximately equal, are measured 
at the time r when the vortex dipole first crosses the line 
y = z = 0. 

Finally, we propose an experiment to probe the physics 
described above. We nucleate vortex-antivortex dipoles 
in the wake of a moving repulsive laser beam which is 
tuned off-resonance from a suitable atomic transition. 
This creates a repulsive potential 
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for the condensate atoms. Here we use r^ — x"^ + [y — 
j/o + vtY , a{z) — CTo-y/l -t- [z/ zjiY and the beam is fo- 
cused to (To = zji/lQ with Rayleigh range z^ = 22 /im 
and Vb = 12.5[1 + tanh(t — T)\hhj. This laser spoon trav- 
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FIG. 2. Crow instability of an elongated vortex dipole in a harmonically trapped Bose-Einstein condensate for d = 285o- (a) 
Kelvin waves are being generated at the ends of the dipole. (b) Two vortex loops have been pinched off, one from each end 
of the dipole. (c) Five daughter vortex loops have been produced through reconnection events, (d) The five vortex loops have 
expanded in diameter appearing as ten vortex lines. The vorticity axis is now aligned along the x axis, (e) A snapshot from 
a simulation using d = 5^o. The x x y x z dimensions of each rectangular box are lOao x lOao x 60ao and the times of the 
snapshots are marked in the frames. 



els in the eondensate for t = 6 ms at a speed v — 0.68 
mm/s corresponding to a Mach number 0.4 and is then 
smoothly withdrawn by ramping down the laser power. 
In Fig. 3a the laser beam is highlighted in the condensate 
density isosurface and its motion is along the y axis. Af- 
ter turning off the laser beam a vortex dipole is revealed, 
which nucleated in the wake of the laser, see Fig. 3(b). 
This vortex dipole becomes spontaneously short circuited 
at its ends due to the combination of the inhomogeneous 
condensate density and the curvature in the laser inten- 
sity profile. The mother vortex loop breaks into five pri- 
mary daughter loops, shown in Fig. 3(c) and Fig. 3(d), 
which undergo multiple reconnections before being con- 
verted into sound waves, which become visible as ripples 
on the condensate surface as time progresses, see frames 
(e)-(g) in Fig. 3. Increasing the Mach number of the trav- 
eling laser spoon increases the number of mother vortex 
rings produced in the wake of the laser. The large number 
of condensate and laser parameters readily accessible in 
experiments enables controlled production of elongated 
vortex dipoles and daughter vortex loops. 

Vortex dipoles and their dynamics in Bose-Einstein 
condensates have recently attracted considerable interest 
[531 132 1311 EZl I3H] • Experiments have considered effec- 
tively two-dimensional systems where the vortex dipoles 
are long-lived structures. We have shown that in elon- 
gated harmonically trapped condensates where axial vor- 
tex degrees of freedom are active, vortex dipoles become 
susceptible to the Crow instability mechanism. Self- 
generated Kelvin-wave excitations cause vortex dipoles 
to disintegrate forming a sequence of vortex loops which 



eventually decay into sound waves. We have proposed 
and simulated an experiment to study this phenomenon. 
Due to the high degree of controllability in this vortex- 
dipole creation method, it could also be used to repro- 
ducibly generate loopy vortex states and to study their 
decay. 

Our results have further implications for decay 
mechanisms of quantum turbulence in harmonically 
trapped Bose-Einstein condensates. In uniform systems, 
a reconnection-dominated Kolmogorov-Richardson cas- 
cade has been suggested to be transformed to a Kelvin- 
wave cascade when vortex-vortex collisions become too 
infrequent to yield a sufficient vortex reconnection rate. 
However, the self-generated vortex loop instabilities may 
be able to sustain vortex reconnections down to the dis- 
sipation scale. It will be particularly interesting to apply 
vortex-dipole nucleation methods to spinor BECs where 
the evolution of non-Abelian vortex dipoles with frac- 
tional charge may lead to drastically different vortex dy- 
namics. 

I thank David Paganin for many useful discussions and 
comments on this manuscript. 
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FIG. 3. Vortex dipole creation and instability in a harmonically trapped Bose-Einstein condensate, (a) A laser beam shown 
in green translates in the y direction, (b) Vortex dipole has formed in the wake of the laser spoon, (c) Decay of the dipole 
has begun, (d) Five daughter loops have been created, (e) Daughter vortex loops have been relinked in further reconnections. 
(f) A snapshot with two vortex rings, (g) Initial vortex dipole has disintegrated creating sound waves seen as ripples on the 
condensate surface. The x x y x z dimensions of each rectangular box are lOao x lOao x 60ao and the times of the snapshots 
are marked in the frames. 
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